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Abstract 

A general method for obtaining the oscillation periods of the interlayer ex- 
change coupling is presented. It is shown that it is possible for the coupling 
to oscillate with additional periods beyond the ones predicted by the RKKY 
theory. The relation between the oscillation periods and the spacer Fermi 
surface is clarified, showing that non-RKKY periods do not bear a direct 
correspondence with the Fermi surface. The interesting case of a FCC(llO) 
structure is investigated, unmistakably proving the existence and relevance 
of non-RKKY oscillations. The general conditions for the occurrence of non- 
RKKY oscillations are also presented. 
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I. INTRODUCTION 



Oscillatory exchange coupling between metallic magnetic layers across a non-magnetic 
spacer has been intensively studied over the last few years. One of the main issues in this 
field has been the determination of the oscillation periods of the coupling as a function of the 
spacer thickness. The physical mechanisms proposed for explaining the so-called oscillatory 
interlayer coupling include the quantum well theory (QWT) of Edwards et al. |I],0|§ an d an 
extension of the Ruderman-Kittel-Kasuya-Yosida (RKKY) theory adapted to the multilayer 
geometry due to Bruno and Chappert 0,^]. 

Both theoretical approaches seem to agree that the origin of the oscillation periods is 
in the structure of the spacer Fermi surface (FS). In fact, in the RKKY theory the periods 
are associated with wave vectors directly obtained from the FS. More specifically, they are 
given by vectors perpendicular to the layer planes which span two points of the FS with anti- 
parallel Fermi velocities |4]|| . There has been a general belief that all oscillation periods are 
given by the RKKY theory and therefore that such a simple picture is always valid. Indeed, 
it has been shown that in certain simple models the periods predicted by the QWT 

coincide with the RKKY ones and are given by the extremal dimensions of the spacer FS 
in the direction perpendicular to the layers. In another case ||10|| , where the lattice lacks 
reflection symmetry about the layer planes, the correspondence between the quantum well 
and the RKKY periods is rather subtle but still obtains. In the models mentioned above, 
harmonics of the RKKY periods are obtained but no new fundamental periods. Furthermore, 
d'Albuquerque e Castro et al. || showed analytically for a very general model that RKKY 
theory holds in the limit of small exchange splitting in the ferromagnetic material. 

However, using a more general approach based on the existence of quantum well 

states |],|2],|3[] , fundamental periods not predicted by RKKY theory have been recently ob- 
tained for the particular case of a FCC lattice along the (110) direction In that commu- 
nication, within the single-band tight-binding model with nearest neighbour interaction, it is 
shown that in addition to the RKKY predictions, a number of distinct periods is involved in 
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the coupling. Those new periods arise from non-negligible higher order terms which RKKY, 
as a second order perturbative approach, does not account for. The resemblance of the FS 
considered with those of the noble metals suggests possible real effects in systems having 
these materials as spacers. Moreover, the occurrence of non-RKKY behaviour within the 
simplest single-band model indicates that it should also happen in the more elaborate multi- 
band case. In fact, non- magnetic transition metals with their multi-sheet Fermi surfaces are 
suitable candidates for presenting non-RKKY periods. Furthermore, in most cases inves- 
tigated so far, each oscillation period comes from distinct q||-points of the 2-dimensional 
Brillouin zone (BZ). Even when multi-periodic oscillations occur, one value of qy yields only 
one period. The existence of more than one period for a single value of qn, which may 
arise in a number of cases, has not been thoroughly investigated and the combination of 
such periods may lead to new features of the oscillations. For the reasons above, a general 
understanding of the mechanisms determining all oscillation periods is needed. 

The purpose of this work is to generalize the method previously used to calculate the 
oscillation periods of the interlayer coupling. The relation between the periods and the 
spacer FS is clarified and in addition to the determination of the usual periods, the general 
conditions for occurrence of non-RKKY behaviour are also established. Furthermore, here 
we present a rather detailed set of results in order to illustrate better the different stages of 
the calculations. The sequence of the paper is as follows. In section || we present a general 
method for obtaining the oscillation periods of the interlayer coupling, where in addition 
to the determination of the usual periods, the existence of non-RKKY behaviour is proved 
to be possible. Similarly to previous theoretical approaches [^,^|,|I|,p|,[16| , an expression for 
the coupling in the asymptotic limit of large spacer thickness is derived. Here however, the 
occurrence of multi-periodicity for a single value of qii becomes evident at the same time 
that the relation between the oscillation periods and the spacer FS is clarified. In section 
PTI] , in order to illustrate the method, we look at the particular case of the FCC lattice 
along the (110) direction which unmistakably shows the relevance of terms beyond RKKY 
approximation. Finally, in the last section we conclude and present the general conditions 
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for occurrence of non-RKKY oscillations. 



II. GENERAL METHOD 

Based on the formalism introduced by d' Albuquerque e Castro et al. 0, where the 
coupling is expressed in terms of the one-electron Green's functions of the multilayered 
system, we here present a general method for determining all the oscillation periods involved 
in the coupling. 

For the sake of simplicity, we consider a system consisting of two parallel magnetic planes 
embedded in an infinite non-magnetic material. We label the two ferromagnetic planes 
and n, so that the number of non- magnetic spacer planes is n — 1. As far as the coupling 
as a function of the spacer thickness is concerned the number of magnetic planes does not 
affect the periods ||, influencing only phase and amplitude of the oscillations. Since the 
periods are our main concern, the restriction to two single magnetic planes does not pose 
any limitations on the results here obtained. 

The expression for the coupling AQ(tc), namely the difference in the thermodynamical 
potential Q between the ferromagnetic and antiferromagnetic configurations, is given by || 

Afi(vr) = ±£ [duf(u)F(q h u) , (1) 

where the function JF is given by 

^(qil, w) = /mtr ln{ l+4Gyq||,w)y« e Gj n (q||,w)K*}. (2) 

In the equations above G£ (q||, uS) is the propagator between planes and n for electrons with 
spin er in the ferromagnetic configuration of the system, f(u) is the Fermi function, V ex is a 
diagonal matrix in orbital indices representing the strength of the local exchange potential 
in the ferromagnetic layers, and the summation over qn is restricted to the two-dimensional 
Brillouin zone (BZ). It is worth stressing that all matrices above are in orbital indices. We 
assume for simplicity that the non-magnetic part of the potential in the ferromagnetic planes 
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is identical to the spacer on-site potential. This is not an essential assumption and can be 
lifted without major effects on the results. In fact, it has been shown that the presence 
of a spin-independent square well does not influence the periods of the interlayer coupling, 
affecting only the phase and amplitude of the oscillation ||. As such, the exchange splitting 
of the magnetic material acts as two localized perturbations on the planes and n and in 
this situation the propagators G can be written as a function of the bulk spacer propagators 
g$ n and g n0 , i.e., 

G no(q||^) = G ° (W% w )>0no(% w )) • ( 3 ) 

Within the particular case of the single-band model the explicit dependence of G G upon the 

off-diagonal propagators go n and g n o is, for each qy and u, given by 

n U _ i T hi9no(l-T hi g 00 ) 
V ex {1 - r hl g^T^t g 0n ) 

where gg m is a general matrix element of the bulk spacer Green's function and = 
±V^j;(l iVe^goo) 1 - An analogous expression can be obtained for G^- Note that whereas 
Eq.([3[) involves matrices in orbital indices, all terms in Eq.(|j) are scalars. 

Returning to the general case and for reasons of future comparison, it is worth presenting 
the RKKY-limit expression for the bilinear coupling J± [ill- It is given by Eq.(|f|) with 
T replaced by 

Frkky{<\\\,oj) = Imtr { 2 V? x g n0 (q\\, u) g 0n (c[\\, u)j . (5) 

We recall that the expansion of AQ(6) in powers of cos(6 l ), being 9 the angle between the 
magnetization vectors in the magnetic planes, defines the bilinear and biquadratic couplings 
Ji and J2, respectively Moreover, it is found that higher order terms than cos 2 (8) are 
negligible so that J x = AQ(n)/2 0. 

It is clear from Eqs.(0), ([3D and (^) that in either case the spacer thickness dependence 
of the interlayer coupling is entirely determined by the off-diagonal propagators c/on and gno- 
The investigation of the oscillation periods thus requires an analysis of g n o as a function of 
n. 
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A. Calculation of the 1-dimensional propagators 



For fixed values of qn , the Green's functions g are equivalent to 1-dimensional propagators 
with qy-dependent hoppings and in-plane energies. We recall that g is a matrix over orbital 
indices and therefore a general matrix element represents the propagation of an electron 
with orbital \x from plane I into an orbital v at plane m. It is given by 



where q± is the wave vector perpendicular to the layers, d is the interplane distance, 
E s ((l\\, q±) describes the bulk spacer band structure, being s the band index, and a s ^(q±) = 
(q± s \q± /i) is the projection of the eigenstate s into the orbital \i for a given q±. 

The integral in Eq.(||) can be evaluated for i < m by simply extending q^ to the complex 
plane and changing the contour integration from a straight line on the real axis to the 
boundaries of a semi-infinite rectangle in the upper-half plane whose base lies on the real axis 
between —ir/d and 7r/d. For the case i > m, the extension is to a rectangle in the lower-half 
plane. In either case the integrand vanishes as \Imq±\ — > oo and because q± = —^+iy and 
q± = ^+iy are equivalent wave vectors, the integrals along the vertical sides of the rectangle 
cancel each other, simplifying the problem to the calculation of the residues associated with 
the poles of the integrand inside the closed contour. The poles are in turn given by the 
values of q± satisfying the condition u + — E s (q\\, q±) = 0, where u + = uj + i + . 

The poles are clearly dependent on the band structure. Due to the in-plane symmetry 
of the system, the general tight-binding Hamiltonian of the homogeneous system is written 
as 



where if rr /(q||) are operators describing the electron hopping between planes r and r'. The 
upper limit of the sum depends on the maximum number of planes f connected through the 
electron hoppings as well as on the lattice structure. As we shall see, in the FCC lattice 






(7) 
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along the (110) direction for instance, even with only hoppings between nearest atoms being 
included, electrons may hop between nearest and next nearest planes. For a fixed value of 
u + , the wave vectors q± satisfying the equation u + — E s (q\\,q±) = are the solutions of a 
polynomial equation in e iq±d , i.e., 

S T r(qn) eiq±rd = » (8) 

r 

where T r (q||) is a function of the qy -dependent hoppings and in-plane energies. Note that the 
solutions of the equation above are exponentials e tq±d and not wave vectors q±_ themselves. 
We label the poles q?, where j numbers the poles of a given band index s. Since E s (q±) = 
E s {—q±) poles appear in pairs and may be real or complex. We see from a simple Taylor 
expansion of uj— E s (q\\, q±) around the real poles that, for a given pair, which pole contributes 
to the integral depends on the sign of the derivative (dE s /dq±) q s. For a given band s which 
crosses the energy w twice, yielding two pairs of real poles, the contributory poles q{ and gf 
have different signs. In other words, when q{ is in the range [0, ir], gf is in the range [— n, 0] 
or vice versa. 

It is interesting to look at the physical significance of the real poles. Being obtained from 
the band structure E s (q\\,q±), the poles q? are just the wave vectors with which electrons 
of energy uj having in-plane energies e(qy) and hoppings t(qy) propagate across the spacer. 
Thus, for uj = Ep, where Ep is the Fermi energy, the real poles are just the perpendicular 
coordinates of the FS for a fixed value of qy. As discussed later, complex poles contribute 
significantly only in the case of thin spacers. 

By adding the residues associated with the contributory poles we obtain an analytical 
expression for gg m , given by 

Wqii,") = EE^e~ i ^~ m)d , (9) 



where the matrix elements of the matrix A S j are 



A%(q h w) = -ida* s ^)a sv {q^) 



. -i 

d£ s (q||,g±) 
dq± 



(10) 



The coefficients Ag% depend neither on I nor on m and the only dependence on the distance 
between the planes is inside the argument of the exponentials. In addition, the exponentials 
are independent of the orbital indices, which means that all matrix elements oscillate with 
the same periods. Note that in calculating Eq.(^) we have assumed that the coefficients 
a Sfl {q_[_) are not ill-behaved, that is, they have neither singularities nor branch points inside 
the integration contour. In this way, we just need to evaluate the coefficients at the poles 

Nevertheless, in a general case the coefficients a Sfl (q±) are not well behaved. Although 
they do not have singularities like in the denominator of Eq.(^j), branch points and cuts do 
exist. The contributions from the branch points must be taken into account otherwise a 
simple summation of the residues will not give the correct result of the integration. This 
is a very important point that, if overlooked, may lead to erroneous results. In fact, in a 
paper by Bruno (17]] where the coupling is expressed in terms of transmission and reflection 
coefficients of electrons across quantum barriers, a similar integral to that in Eq. (|6|) arises. 
In that communication no particular attention was given to the analytic properties of the 
coefficients a Sfl (qj_) and the result above was erroneously concluded to be exact for arbitrary 
values of i and m. It is important to stress that Eqs.(^) and ( |TUD are not correct in general. 
The equations above become exact for all I and m only within the single-band model, where 
there are no coefficients introducing branch points. This is clearly seen in Fig.|l] where the 
real and imaginary parts of the function posO^) °f a linear chain with hoppings up to four 
nearest neighbours within the single-band model are displayed. The parameters used were 
e = 0, t\ = —1/2, t 2 = —1/2, t 3 = 1/5 and i 4 = 1/10, where e is the in-plane energy and 
tt are the hoppings to the £-th neighbour planes. The analytical results following Eq.(§), 
represented by the lines, are in absolute agreement with the brute-force calculation of Eq.(^), 
represented by the symbols. In the general case there are also contributions to the q± integral 
arising from branch cuts associated with the coefficients a SAt (g_i_). Since in general these cuts 
do not intersect the real axis, the factor e l <}±( i - m ) d j n the integrand of Eq.(^) ensures that the 
contribution is negligible for large \£ — m\. Thus, failure to treat the branch cuts correctly 
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does not introduce errors in the asymptotic limit of large \£ — m\. 

It is clear that gi m is an oscillatory function of the distance between the planes and it 
oscillates with different periods, each one associated with its respective wave vector gj. In 
addition, the combination of exponentials allows a quasi-periodic behaviour of ge m because, 
in general, the wave vectors are incommensurate. This fact, as we shall see, may have a 
striking effect on the oscillatory coupling. Another point to be mentioned is the possibility 
of evaluating <% m in Eq.@ for non-integer values of I and m, which allows the evaluation of 
the coupling for continuous values of spacer thickness. 

We now mention the case of complex poles. Since they have non-zero imaginary parts, 
their contribution for large separation between the planes is again negligible. Thus only 
real wave vectors contribute to the oscillations for large \£ — m\. Bearing in mind that the 
principal contribution to the coupling comes from the energy region around the Fermi level 
Jl|fJJ^,|l^j , it is clear that at least one real wave vector will be involved in the expression for 
ge m . Thus the sum over the band index s in Eq.(|9]) must be only over the bands which cross 
the Fermi level. 

B. Multiple Fourier expansion 

Having investigated the function g^ m as a function of the distance between the planes, we 
recall that T is expected to oscillate with the same periods. The obvious way of representing 
the function T is through a simple Fourier expansion, which nevertheless cannot be done 
due to the quasi-periodicity of g. In other words, because the wave vectors gj are in general 
incommensurate, go n and g n Q may have an overall non-periodic behaviour. To overcome this 
problem we make use of the multiple Fourier series, recently used to deal with the problem 
of magnetic thickness dependence of the interlayer coupling || . The procedure consists in 
replacing the spacer thickness n multiplying the wave vectors gj in each exponential with 
fictitious spacer thicknesses n|, such that Eq.(|9|) becomes 
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The new function g is now periodic in each variable separately, and so is its corresponding 
T . Bearing in mind that g does not correspond to the physically real case, which is repro- 
duced by making n| = n for every j and s, we evaluate the Fourier series for each variable 
rij separately, eventually making them all equal n. T then becomes 

T E C^ m ^e ini ^ m ^\ (12) 
where are integers, C^ m i^ m s^ are the Fourier coefficients given by 

C K ,... lW |,...) = (jT^y-) C dn l- dn l- ^K,...,^,...)e id ^K^"5) , (13) 

and = 27r/|g|| are the periods associated with the wave vectors gj. The number of indices 
involved in the Fourier coefficients is the same as the number of poles because each q s - has 
been associated with a fictitious spacer thickness nl In spite of the somewhat congested 
notation, in practice there are only a few indices in the coefficients due to the restricted 
number of real solutions crossing the Fermi level. In fact, within the single-band model 
the calculation of the Fourier coefficients can be extremely simplified, as shown in appendix 
A. The analysis of such a simplified model, although primarily illustrative, reflects some 
properties common to the more general case. Two of these properties are C_ mi _ m2 = C^ i m2 
and C mi>m2 = for odd values of m x + m 2 . 

1. RKKY- Coefficients 

We now focus on the RKKY-limit expression for the coupling. This is the limit of small 
ferromagnetic exchange splitting. It is clear from Eq.([|) that unlike the coupling Af2, the 
RKKY expression involves neither logarithmic functions nor combinations of G^ and Gq u . 
In fact, J 1 rkky simply contains a single product of the propagators g n0 and g 0n . Bearing 
Eq.([5p in mind we define Trkky = [J~'rkky ~ (^rkky)*]/^ an d from the expression for an 
arbitrary element gi m in Eq.(|9|) we write 
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Prkky = 2 K'EEE 4? K*y ^ +q > d • (14) 

s,s' 

Clearly, the above equation shows that the integrand Trkky can be manipulated and put in 
a Fourier expansion form, similarly to T in Eq. (|T2|) . However, the main difference is that in 
this case the number of terms involved in the expansion is limited, related to the maximum 
number of poles in Eq.@, whereas the sum in Eq.(|l^) runs over an infinite number of 
terms. Such difference is due to the fact that the RKKY-limit expression results from a 
perturbative approach corresponding to the lowest order term of an expansion. All the 
terms included in Trkky are also included in the full expression T . The higher order terms 
correspond to oscillations that are either harmonics of the fundamental RKKY behaviour 
or effectively new oscillation periods. This point confirms the RKKY theory as a limit of 
the more general approach used here. We also point out that in this case, the argument of 
the exponentials has the form (mjgj + m^g|/'), involving only two wave vectors and with the 
constraint |m| + m s -,\ = 2, being \m^\ < 2 and \m s -,\ < 2. This is a consequence of Frkky 
being the product of only two propagators, namely, g n o and go n . In the more general case 
of J 7 , the infinite possibilities of combining the mentioned propagators does not pose such a 
restriction. 

At this point we raise the question of whether the RKKY periods of the interlayer 
coupling remain valid as the exchange splitting of the ferromagnetic layers, associated with 
V ex , is increased. The question can be reformulated as whether the higher order terms in 
J 7 , apart from the RKKY harmonics, do contribute to the coupling. The answer to this 
question is in the evaluation of the integrals in Eq. ([[]). 



C. Stationary phase approximation 

Having written T (and Frkky) in terms of a Fourier expansion, we can now make 
use of the stationary phase method [|^p|,|3|,|T6|| > which yields an analytical expression for the 



coupling, asymptotically exact for thick spacers (n » 1). 
Inserting Eq.([L2|) into ([[]), we have 
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A<>U) = - £ E / ^/H^,„.) e ^^K«P . (15) 

(mj,..,mj,...) q|| 

We recall that the coefficients C( m i v .. )m s ...) and the wave vectors gj depend on both integra- 
tion variables qy and uj. The stationary phase method shows that in the asymptotic limit 
of large spacer thickness the only contributions to the integral over qn arise from the neih- 
bourhood of points which make the arguments of the exponentials stationary. We call these 
points qjp where the superscript a labels the different vectors. Furthermore, the integral 
over uj is also non- negligible only around the Fermi energy Ep. Hence, the wave vectors that 
effectively contribute to the coupling must satisfy the following condition, 

E[^V||^(q||,^)] = 0, (16) 

where Vy is the two-dimensional gradient in qy -space. In other words, the equation above 
gives the necessary conditions for a constructive interference between the electrons across the 
spacer. The otherwise destructive interference does not contribute to the coupling. Eq . (|16|) 
simply selects the effective periods with which the coupling oscillates. 

Note that for uj = Ep, each wave vector gf (qy, Ep) represents a surface in the reciprocal 
space and when put together, the surfaces gj reproduce the spacer Fermi surface. Following 
this argument and bearing in mind that non-real values of g| strongly damp the oscillations, 
we stress that only real solutions must be taken into account. The possible stationary 
solutions of Eq.(|T6"D reflect the lattice structure and are usually located at points of high 
symmetry, even though they may occur at general points of the 2-dimensional BZ. 

In such limit we obtain an analytical expression for the coupling given by [1],0,§ 

A^)=E E % „r <''"'"• (17) 

a (mj,...,m^,...) 



fc-(m{,...,m°,...) y _ ));/ J 



/ I- ,\ rr a r a I d 2 g d 2 q \-l/2 

I 2\/2d\ a (ml,...,m s j ,...)^(m 1 1 ,...,m s j ,...)\ dql dq^ 



sinh 



(18) 



where q = J^sji^Qj ): Qj a = QjiQp Ep), and (3 = 1/ksT. Here we neglect the energy 
dependence of the Fourier coefficients which in some cases |16] can be important. The factor 
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a? i s \ takes the value i for a minimum, — i for a maximum and 1 for a saddle point of 

(m^,...,rnj,...) ' " 

the surface q(q x , q y ). The sum over a regards all different stationary solutions and the inner 
sum over (m\, m|, ...) runs over the values allowed by the stationary phase condition. The 
coupling oscillates with periods X? i « v = 2tt/\ J2si( mS iQi a )\ an d their weights depend 
on geometric factors as well as on the Fourier coefficients. 

An important issue is the relationship between the Fermi surface and the wave vectors 
that effectively contribute to the coupling. The lowest order combinations amongst the 
possible wave vectors ?!(q|[, Ep) bear a direct correspondence with the FS. Indeed, this cor- 
respondence agrees with the RKKY geometrical picture for selecting the oscillation periods 
of the interlayer coupling. However, in the general case of a higher order combination of wave 
vectors the RRKY relation no longer holds, and the effective wave vectors are not directly 
related to the spacer FS. This is a very important result which shows that the geometrical 
criteria used for selecting the oscillation periods of the interlayer coupling are incomplete. 
Even though the wave vectors g| are directly related to the FS, the general wave vectors 
J2 S j{ m j1j) are n °t- Furthermore, note that if g„o(qfj\ Ep) has only one oscillation period, 
the full expansion in Eq. ( |TTD contains merely harmonics of the fundamental RKKY periods. 
In other words, if the FS is cut in a single point by the line drawn from the origin perpen- 
dicular to the layers for a fixed qy, the RKKY periods and its harmonics are the only ones 
present in the coupling. In fact, it explains the agreement between the periods predicted by 
the RKKY and QWT in the models previously mentioned |Hflj3|,|10[1 ■ For a multi-periodic 
g n o(q\\, E F ) though, where the FS is intersected in more than one point corresponding to 
incommensurate wave vectors, the higher order combinations of the wave vectors correspond 
neither to a RKKY period nor to any of its harmonics, but to fundamentally new periods. 

Multi-periodicity of the off-diagonal propagator is a necessary but not sufficient condition. 
The argument above simply shows that non-RKKY behaviour is possible in presence of 
multi-periodic g n o. To find out whether it really occurs one has to solve Eq. ([16|) and check 
whether there are real wave vectors q? for values of m\,...,m S j, ... other than the ones obtained 
by the RKKY theory. 
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In summary, the determination of the oscillation periods results from the analysis of 
the bulk spacer off-diagonal propagators where the relevant real wave vectors, under the 
stationary phase condition, are obtained from the band structure. The relation between the 
spacer FS and the oscillation periods is the same as the geometrical RKKY picture for the 
lowest order terms. Nevertheless, for higher order contributions a direct correspondence with 
the FS is not evident and the stationary phase condition displayed in Eq.(|l6|) determines 
the periods with which the coupling oscillates. We have shown that oscillation periods other 
than the ones predicted by the RKKY theory may arise. To illustrate the method and 
unambiguously prove the existence and relevance of periods other than RKKY ones, we 
next evaluate the particular case of the FCC lattice along the (110) direction within the 
single-band tight-binding model. 



III. APPLYING THE METHOD 

Having presented the general method for determining the periods of the interlayer cou- 
pling, here we apply it to the particular case of a FCC(llO) system within the single-band 
nearest-neighbour tight-binding model, which unmistakably shows the existence of non- 
RKKY periods. Moreover, in the light of the QWT this system is interesting because of the 
difficulty in determining analytically the energies of the resonances and size quantized states 
in such a structure ||12|| . For a single band the notation introduced in the previous section 
becomes simpler due to the omission of the band index s. In addition, all the matrices over 
orbital indices become scalars. 

We label the nearest neighbour hopping t an d with no loss of generality assume the 
on-site potential eo = 0. The bulk density of states of the FCC lattice is displayed in Fig.||. 
It is worth pointing out that there are different contributions to the interlayer coupling, 
depending on the position of the Fermi energy inside the band. Three distinct regions are 
observed and their boundaries are shown by the vertical dashed lines in Fig.[| In addition 
to the traditional polyhedral BZ, we make use of an auxiliary prismatic BZ due to the 
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2-dimensional symmetry of the multilayers ||. Both zones are equivalent and enclose the 
same volume in reciprocal space. The energy boundaries indicated correspond to Fermi 
energies such that the FS just touches the two distinct Brillouin zones. More specifically, 
for Ep = —4to the corresponding FS touches the prismatic BZ, whereas Ep = is the exact 
Fermi energy where the FS opens the necks which are familiar from the noble metals W\ . 



As previously mentioned, even within the nearest neighbour tight-binding model, in this 
case electrons are capable of hopping between nearest and next nearest planes. The bulk 
band structure is then given by 

S(q||,?j.) = e(q||)+2ti(q||)cos(gj.d)+2t 2 (q||)cos(2gj.d) . (19) 

We assume the z direction perpendicular to the layers. The 2-dimensional BZ over which 
the stationary solutions must be searched for is a rectangle defined by —7r/2d < q x < ir/2d 
and ~Tc\^2/Ad < q y < ny/2/Ad, where q x and q y are the components of qn. The explicit 
dependence of E on qy is given by e(qy) = — 2 1 cos(2q x d), ti(q||) = — 4t cos(q x d) cos(y/2q y d) 
and £ 2 (q||) = —to- 

Following the steps of the previous section, u + — E(q\\,q±) = gives the poles which 
yield the periods of g n Q. Eq.([T9|) can be rewritten in terms of a second degree polynomium 
in cos(g_i_<i) whose solutions qi and q 2 are given by 



cos(g i d) = - ( 2 - \ , (20) 

where 7(qy) = t 1 /t 2 and j is either 1 or 2. Nevertheless, the poles qi and q 2 are in the 
g^-complex plane. By mapping the relationship between the two complex functions we can 
uniquely relate one value of q± inside the integration contour to the corresponding value 
of cos(g_i_<i) and vice versa, as shown in Fig.^. The equation above shows that cos(gi d) 
and cos(g2 d) have opposite imaginary parts and therefore, following the maps in Fig.|3], qi 
and q 2 have opposite real parts. This confirms what has been said in the previous section 
concerning the position of two different wave vectors for a given band s. 

In addition, looking at Eq.(f2l]) with uj = Ep it is clear that non-RKKY behaviour cannot 
arise in the bottom region of the band (— 12to < Ep < —4to). This is due to the fact that 
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in such a region q± and q 2 cannot both be real. In fact, the corresponding FS in this region 
is spherical-like. However, in the range — 4t < Ep < 4t , corresponding to the top regions 
inside the band, both q\ and q2 may be real indicating that non-RKKY periods are possible. 
Having obtained the poles, from Eqs.([5|) and ([H]) we have 

<7 n0 (q||,a,) = A 1 (q h u) e ^H>"W + A 2 (q h to) e^^ nd , (21) 

where 

Aj(q\\,u) = ^2i [cos(gi<i) — cos(g 2 ^)] \Jl — cos 2 (g.,<i)j . (22) 

In the present case go n = g n Q. Eq.(|2T|) shows that for values of qy and u, for which qi and q 2 
are real, g n o oscillates with the superposition of two periods, 27rd/\qi\ and 2ird/\q2\, which 
are in general incommensurate. In those cases g n o exhibits a quasi-periodic dependence on 
n. 

Because there are only two poles in this case, the function J 7 , following Eq.([L2|), becomes 



JT — Q e i(miqi+m2q2)nd (23) 

rai,m2 



where the Fourier coefficients are given by Eq. (|TB|) . Within the single-band model, the 
Fourier coefficients can be calculated in an alternative way, shown in Appendix A. Following 
the Appendix, we recall that the lowest order Fourier coefficients correspond to the RKKY- 
limit J-rkky- They are C mitm2 where m!+m 2 = 2, \mi\ < 2 and |m 2 | < 2. More specifically, 
C 2 ,o = C ,2 = C 1A = and the property C_ mi _ m2 = (C mi , m2 )* define all the 

six coefficients present in the RKKY expression, where T = 4r|Tj (1 — (7oo)(1 — r i 9oo) an d 
T|| are defined after Eq.(|). Beyond the lowest order coefficients are the terms expected to 
contribute to the non-RKKY behaviour. It is then illustrative to present one of these terms, 
namely C 3 u 



W + (nr - ¥ 



(24) 



The stationary phase condition in Eq.(|T^) is rewritten as 
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miV|| ft(q||, E F ) + m 2 V|| q^(q h E F ) = . (25) 

The condition above imposes that both gradients must be colinear and moreover, the ratio 
between their moduli must be a rational number. The determination of the stationary 
points qjj* = (q x ,q y ) yielding the periods A^ im2 = 2n/\m 1 q^ + m 2 q^\ is rather tedious 
but straightforward. As previously mentioned, the solutions are often located at highly 
symmetric points of the BZ but more general solutions may occur. Nevertheless, because 
the RKKY periods as well as the most important higher order terms come only from those 
high-symmetry points, we shall firstly focus on these solutions. Across the whole energy 
band, four stationary points are found on the corners of the 2-dimensional irreducible BZ. 
They are qj, = (0,0), qj = (0,±7iV2/4gT), qjf = (±7r/2d, =br-\/2/4d) and qjf = (±vr/2d,0). 

In the first region A^q is the only period, which is actually present along the entire 
band. The graphical relation between the spacer FS and the oscillation period is displayed 
in Fig.[|, where the FS for an arbitrary energy (E F /2t = —3.0) is shown together with the 
wave vector g" and the cross section of the prismatic BZ for q y = 0. The cross section is 
defined by —7i/2d < q x < ir/2d and —ir/d < q z < ir/d. Note that 2q" corresponds to the 
extremal wave vector spanning the FS, in according to the RKKY picture. Clearly, higher 
order periods A TOl) o for mi > 2 are merely harmonics of the fundamental A2,o- 

In the second region an additional period \\ 1 arises. Similarly, the vector q\ + q\ yields 
the fundamental period and possible higher combinations mi(q\ + q\) add no new terms but 
harmonics. Fig.[| shows a cross section of the FS (E F /2t = -1.0) for q y = ny/2/Ad. The 
wave vectors qj and q^ are highlighted and it is clear that qj+q^ satisfies the RKKY criteria. 
At first glance one could think, accordingly to Fig|| that 2qJ and 2q^ are also contributory 
wave vectors for the same reasons as in the first region. Nevertheless, by looking at the 
g y -dependence of the FS we find that those vectors are not solutions of Eq. (|25|) . 

The most interesting case occurs in the top region, where the FS has necks resembling 
those of the noble metals. In addition to A^q at qff, other stationary solutions occur at qjj 
and q[i . At qf?, analogously to the second region, the solution comes from qf + q$ and does 



17 



not introduce new periods. At qjj though, both gradients in Eq.(p5[) vanish simultaneously 
allowing any values of mi and m 2 to satisfy the stationary phase condition. Therefore, the 
coupling oscillates with as many periods m2 as there are non-zero coefficients C^ m2 . 
The lowest order A 20 , Aq 2 an d XI i correspond to the RKKY periods and the geometrical 
relation between the wave vectors q£, q 2 and the FS in Fig.|6| confirms this point. For the 
present model we find that q{ — q^ = n. Thus, the oscillation periods A 2 and Aq 2 cannot be 
distinguished just by looking at discrete integer values of n. In addition, AJj = Xf 1 , which 
leads the third energy region to contain three non-equivalent RKKY periods. We recall that 
in FCC Cu Ep lies in the third energy region with FS necks. A quantitative description 
of the oscillation periods for Cu can be obtained within the present framework by going 



beyond nearest neighbours and using the tight-binding parameters of Halse ||13|| . Then A \ x 
and Xf x become distinct periods and, by taking into account the equivalence of A 2 and Aq 2 
for a discrete lattice, we find exactly the four RKKY periods of Bruno and Chappert Q. 

However, because in general the wave vectors qj and q 2 are incommensurate, higher order 
periods are not simple harmonics of the RKKY ones, unlike in the two lower regions. In fact, 
the period A3 x does not correspond to any harmonics of the fundamental RKKY periods. 
Furthermore, even though the wave vectors qj and q 2 separately have clear correspondences 
with the FS, as shown in Fig.|], the wave vector 3qJ + q 2 has not. From this argument it 
becomes clear that higher order periods do not satisfy the geometrical picture of the RKKY 
theory, having instead an indirect relation between the wave vectors and the FS. 

Fig.^ summarizes the results obtained here by displaying the oscillation periods as a 
function of the Fermi energy Ep. The full lines labeled from 1 to 5 are the fundamental 
RKKY periods A 20 , Aj v A 20 , Aq 2 and AJ 1; respectively, and the dot-dashed line in the top 
region represents only one of the non-RKKY periods, namely A^. Recalling the relation 
Qi ~ <?2 = we see t na t the oscillation periods A 2 and Aq 2 cannot be distinguished just by 
looking at discrete integer values of n. Still focusing on the top energy region, Fig.[S] exhibits 
the bilinear coupling J\ = AQ(tt)/2 at temperature T = 2.0 x 10~ 3 W/kB as a function of the 
spacer thickness for E F /2t = 1.64 and V ex = 0.15W, where W = 16t is the spacer band- 



width. This value of Ep is chosen so that A3 1 is well separated from the RKKY periods, 
as may be seen in FigfT]. The full line in Fig.|8| corresponds to the result obtained from 
Eqs.([l|) and (0), whereas the dashed line is the RKKY approximation, following Eqs.(fTJ) 
and (|5|), scaled down by a factor of 8 for reasons of comparison. Despite the agreement 
in the dominant long period, the difference in the fine structure of the oscillations reflects 
contributions beyond the fundamental RKKY periods. 

The evaluation of the stationary phase approximation in this case unmistakably proves 
the existence and relevance of the non-RKKY periods. Fig.|| displays the absolute value 
of the ratio between some Fourier coefficients (7° m2 and C{ x , the largest coefficient of a 
fundamental RKKY period, as a function of V ex . The full line corresponds to the Fourier 
coefficient associated with the new period A3 1 whereas the dashed and dot-dashed lines 
correspond to A^ 2 and A^q, respectively, which are simple harmonics of the fundamental 
RKKY periods. As expected for very small exchange splittings, the magnitudes of higher 
order coefficients relative to that of the fundamental RKKY ones are negligible || , indicating 
that non-RKKY behaviour becomes more important as the exchange splitting is increased. 
An interesting consequence is that, even though the periods are determined from the spacer 
FS, their relevances depend on the nature of the magnetic materials involved. Note that 
for V ex ~ 0.15W 7 the ratio | C3 1 1 / 1 Ci 1 1 ~ 1/2. This indicates that for such a value of 
V ex , which is exactly the same as in Fig. the Fourier coefficients are of the same order 
of magnitude; the only reason why both periods do not have comparable weights is due 
to geometrical factors determined by the curvatures of the respective surfaces. In fact, 
the curvature of the surface qi + q 2 vanishes at qfr, making it impossible to evaluate the 
stationary phase approximation for the period AJ x even though it undoubtedly indicates 
a strong contribution. Taking into account the geometrical curvatures when possible, we 
evaluate the stationary phase approximation and show in Fig.[H] the separate contribution 
of some periods to the exchange coupling. The full line is the new period A3 l whereas the 
dashed line and the the dot-dashed line represent X^q and A^g, respectively. Clearly, A| x is 
as important as the RKKY periods, except for the long dominant one. It is worth pointing 
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out that the rate of decay of the oscillation amplitude associated with the new periods is 
exactly the same as for the lowest order ones, falling off as 1/n 2 for T = 0. Concerning the 
actual amplitude of the oscillations, two points should be highlighted. The strength of all 
periods depends rather critically on the matching of the bands in real materials, which is 
not well reproduced by our single-orbital model, and the omission of the energy dependence 



of the Fourier coefficients could be serious for Cu-based systems [|16 



The solutions treated above all come from the high-symmetry points qy , qjj, qy and 
qj|. Note that fundamentally new periods arise but always in the presence of RKKY ones. 
They result from crossed mixtures of the fundamental wave vectors, which in turn yield 
the RKKY periods. Nevertheless in general, as previously stated, solutions of Eq.(p5[) away 
from high-symmetry points are found and we point out that their contributions are not 
necessarily associated with the existence of RKKY periods. In fact, it can be shown in 
the second energy region that energy- dependent solutions exist along the line q x = with 
corresponding wave vectors qi and q2 which do not yield any possible RKKY period. This 
indicates there are colinear gradients along that line satisfying Eq.Q25|) for different values 
of mi and m 2 even though no RKKY stationary solutions are obtained. These solutions 
are found to be far less significant in their contribution to the exchange coupling than the 
highly symmetric ones. 



IV. CONCLUSIONS 

In summary, we have shown that the RKKY theory is not capable of predicting all 
the oscillation periods of the interlayer coupling, as generally believed. The actual relation 
between the periods and the spacer FS is more subtle than the RKKY geometrical picture. 
Even though the lowest order periods indeed correspond to the RKKY predictions, non- 
negligible higher order terms arise and do not bear any direct correspondence to the FS. 
Instead, they correspond to indirect combinations of fundamental wave vectors, which in 
turn, may or may not yield RKKY periods. The appearance of non-RKKY oscillations 
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clearly requires the multi-periodicity of the bulk spacer off-diagonal propagator as a function 
of the spacer thickness, although this is not a sufficient condition. In addition, real solutions 
of the stationary phase condition must exist in order to have a constructive interference 
between the electrons across the system. The otherwise destructive interference damps the 
oscillation not contributing to the coupling. Finally, if such behaviour can be found even 
within the single-band model, for which the FS is simple and has a single sheet, we may 
expect the occurrence of non-RKKY periods in systems with more elaborate FS. Surfaces 
with more than one sheet, such as non-magnetic transition metals for instance, are possible 
candidates for presenting such periods. In those cases, the interpretation of the results in 
terms of the RKKY theory may be misleading. 

We are grateful to EPSRC and Royal Society of UK, and CNPq of Brazil for financial 
support. We also would like to acknowledge useful discussions with Dr. R. B. Muniz and 
Dr. A. Umerski. 



APPENDIX A 



Alternatively, within the single-band model, the Fourier coefficients C mijm2 in section [El I 
can be calculated in the following way. Bearing in mind the equivalence between g 0n = g n0 , 
Eq.(f|) can be concisely written as 

Gi = r ' 90n %-^ 900) ±(T l9n0 f> , (Al) 

Vex j=o 



and analogously, 



G n = 77 2^ ( T l 9n0) • (A2J 

ex j'=o 



The product of Eqs.(^) and (|A2|) is 

r T #00 ) 



Gi Gl = _ nn(i-^oo)(i-n^) E(TT)2 ; (n) ^' (, n0 ) W+i) . (A3 ) 



Writing the logarithmic function as a series, we have from Eq.(^j) 
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1=1 



.):) 



(A4) 



where T = — 4t|T; (1 — T| g o)(l — r .| Poo)- We recall that the summations above have indices 
with different limits. Whereas the external summation index I runs over the integers starting 
from i = 1 onwards, the inner indices j and j' have zero as their initial values. Eq. can 
be manipulated and rewritten as 

00 



^ = M E ~ — ^ — ^ E - E(n) 20l+ - + ^ ) (n) 2W+ - +j;) (^"••••^•^••••^ 



(A5) 



For a fixed £ there are twice as many inner summations, running from j\ and j[ as far as 
je and j^. We define the integer p as the exponent of the function g^ in Eq . ( |A5|) , i.e., 



p = ji + j[ + ... + ji + j[. Note that p runs over the positive integers greater than zero. 



Squaring g n0 in Eq.(21) and inserting it into Eq. (|^5|) we obtain the function T in terms 
of a Fourier expansion. Clearly, the exponential e n m i , ?i+ m 2<?2)nrf arises from terms of (g n o) 2p 
such that p = mi + m2 , The number of terms in turn depends on the number of possible com- 
binations between the implicit indices j, j' and £, and the larger p is, the more combinations 
are possible. 

For p = 1 for instance, the only possible combination of the indices is j\ — j[ = and 
1=1. The terms coming from this contribution correspond to the RKKY-limit expression 
in Eq.(||), which indicates J-rkky as a subset of the full expression T and confirms the 
RKKY theory as a limit of the more general approach used here 0. 

For p = 2, more combinations are possible. In addition to i = 2 and ji = j[ = = 
j' 2 = 0, four other combinations come from 1=1 and ji + j[ + j% + j' 2 = 1. Note that 
the lowest order coefficients are easily calculated through this alternative method but as the 
exponent p is increased, more terms arise and increase the difficulty. Nevertheless, because 
we are searching for periods other than the RKKY ones, the lowest order terms beyond the 
RKKY-limit are expected to be the important corrections to the coupling. 
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Finally, from the discussion above, two properties of the Fourier coefficients are straight- 
forwardly obtained. They are C_ mi) _ m2 = C* m2 and C mijm2 = unless m 1 + m 2 is even. 
These properties are also valid in the general case of section y. 
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FIGURES 




lm co 



FIG. 1. Comparison between the analytical expression and the brute-force calculation of the 
propagator gos(uj). The lines represent the analytical results and the symbols the numerical ones. 
The full line and the circles are the real part of whereas the dashed line and the squares are 
the imaginary part. The parameters of this tight-binding Hamiltonian are indicated in the text. 
In order to improve the convergence of the brute-force calculations we have varied uj along the 
imaginary axis. 
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FIG. 3. Correspondence between the q± and cos(gj_)-complex planes for £ < m. By considering 
qx_ inside the integration contour (bold path on the left), each point on one plane has an unique 
correspondence on the other one. For the case i > m, since cos(gj_) is even, the map of the 
gx-complex plane must be reflected about the axes. 
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FIG. 4. Cross section of the Fermi surface for Ep = —3.0 (Region 1) and q y = 0. The rectangle 
is a cross section of the prismatic Brillouin zone and the wave vector qf is also displayed. 
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FIG. 5. Cross section of the Fermi surface for Ef = —1.0 (Region 2) and q y = TTy/2/Ad. The 
rectangle is a cross section of the prismatic Brillouin zone. The wave vectors q\ and q\ are also 
displayed. 
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FIG. 6. Cross section of the Fermi surface for Ef = 1.0 (Region 3) and q y = The 
rectangle is a cross section of the prismatic Brillouin zone. The wave vectors qf and q% are also 
displayed. 
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FIG. 7. Oscillation periods of the interlayer coupling as a function of the Fermi energy Ep. 
The vertical dashed lines separate the three distinct regions inside the band. The full lines labeled 
from 1 to 5 are the RKKY periods (see text) whereas the dot-dashed line in the top region is one 
of the non-RKKY periods, namely A| x . 
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FIG. 8. Full line shows the interlayer coupling for Ep/2to = 1.64, V ex = 0.15W and 
ksT = 2.0 x 10 _3 P^, whereas the dashed line is the RKKY prediction for the same parameters 
scaled down by a factor of 8 for comparison. 
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FIG. 10. Separate contributions to the interlayer coupling through the stationary phase ap- 
proximation. The full line is A^, dashed line is A^o an d dot-dashed line is A^o- Ep/2to = 1.64, 
V ex = 0.15W and k b T = 2.0 x IQ^W . 
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